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We study high efficiency, multi-terawatt peak power, few angstrom wavelength. X-ray Free Elec¬ 
tron Lasers (X-ray FELs). To obtain these characteristics we consider an optimized undulator 
design: superconducting, helical, with short period and built-in strong focusing. This design re¬ 
duces the length of the breaks between modules, decreasing diffraction effects, and allows using a 
stronger transverse electron focusing. Both effects reduce the gain length and the overall undulator 
length. The peak power and efficiency depend on the transverse electron beam distribution and 
on time dependent effects, like synchrotron sideband growth. The last effect is identified as the 
main cause for reduction of electron beam microbunching and FEL peak power. We show that the 
optimal functional form for the undulator magnetic field tapering profile, yielding the maximum 
output power, depends significantly on these effects. The output power achieved when neglecting 
time dependent effects for an LCLS-like X-ray FEL with a 100 m long tapered undulator is 7.3 
TW, a 14 % electron beam energy extraction efficiency. When these effects are included the highest 
peak power is achieved reducing the tapering rate, thus minimizing the reduction in electron micro- 
bunching due to synchrotron sideband growth. The maximum efficiency obtained for this case is 9 
%, corresponding to 4.7 TW peak radiation power. Possible methods to suppress the synchrotron 
sidebands, and further enhance the FEL peak power, up to about 6 TW by increasing the seed 
power, are discussed. 

PACS numbers: 41.60.Cr, 41.60.Ap 


I. INTRODUCTION 

Maximising the extraction efficiency of a Free Electron 
Laser (FEL) via undulator tapering has been examined 
theoretically numerically El a 0 and demon¬ 

strated experimentally from the microwave 0 down to 
hard X-ray wavelengths [7]. At hard X-ray wavelengths 
the achieved efficiencies are a factor of 3 larger than the 
exponential saturation value [7]. Single molecule imag¬ 
ing applications require TW-level X-ray pulse power and 
thus an improvement in the extraction efficiency of more 
than an order of magnitude in the next generation of X- 
ray FELs [5]. This has sparked renewed interest in the 
community with a number of recent studies specifically 
devoted to finding the optimal tapering law using a model 
based approach and a form of parametric optimization [1] 
0 0. With the exception of Ref. [S], which considers 
longer wavelengths than we do in this paper, the opti¬ 
mizations presented have dealt exclusively with the time 
independent physics of tapered FELs, producing a taper 
profile which maximises the extraction efficiency ignor¬ 
ing time dependent effects like the synchrotron sideband 
instability [T]. The physical difference between time in¬ 
dependent and time dependent optimization arises from 
noise in the electron beam current distribution and slip¬ 
page of the radiation field. This can drive the amplifica¬ 
tion of parasitic frequencies and the sideband instability, 
causing temporal fluctuations in the electric field profile, 
particle detrapping and eventually taper saturation. 

In this paper we show that the solution obtained for 
the optimal taper profile in time independent simulations 
does not yield the maximum extraction efficiency when 


fully time dependent physics is included in the dynamics 
of the the electron beam and radiation field system. We 
study the optimization problem by following the multidi¬ 
mensional scan method of Ref. [5] for a superconducting, 
2 cm period helical undulator with built in focusing. This 
undulator design is optimized for maximum efficiency, re¬ 
duction of intra module undulator length, strong trans¬ 
verse focusing, short gain length and minimum total un¬ 
dulator length. The characteristics are given in Table 
1. The simulations are performed for a transversely flat 
electron beam distribution as this maximises the output 
power, as described in Ref. nni. 

The paper is divided as follows. In section II we discuss 
the undulator design, outlining the undulator parame¬ 
ters and the feasibility of practical realization. In section 
III we discuss the self-seeding method and simulate the 
SASE FEL upstream of the optimized tapered undulator. 
In sections IVa and IVb we present the time independent 
and time dependent tapering optimization results and 
discuss the differences between the two. In section V we 
analyze in detail the sideband instability, its onset, am¬ 
plification and its impact on particle dynamics and taper 
saturation. Finally we discuss methods to suppress the 
instability and outline future work on high efficiency ta¬ 
pered undulators and optimization algorithms. 


II. UNDULATOR DESIGN 

We apply the tapering optimization method [3] to an 
undulator designed specifically to achiveve TW power X- 
ray pulses in the shortest possible undulator length. Our 
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FIG. 1: Schematic of the undulator for hard X-ray multi TW 
peak power output, designed to achieve high extraction effi¬ 
ciency in the shortest possible distance. 

ideal undulator is superconducting, with a short 2 cm 
period and a peak on axis field Bq of 1.6 T. For a double 
helix bifilar magnet with equal and opposite currents this 
field is given by m- 

Ah T 

Bo = [kuaKo{kua) + iFi(/c„a)] (1) 

where I is the current in the coils, = 2TTjXu is the 
undulator wavenumber, a is the helix radius and Kq and 
Ki are modified Bessel functions. For a helical bore ra¬ 
dius 0=7.5 mm the total current required through the 
coils is J = 484 A which, considering coils of ~ mm^ 
surface area, gives a current density below the critical 
value for superconducting NbTi or NbSSn wires. From 
the point of view of operation a superconducting undula¬ 
tor has advantages such as resistance to radiation damage 
and reduced sensitivity to wakefields, for a more detailed 
description see Ref. [T^]. The undulator is helically po¬ 
larized as this increases the effect of refractive guiding in 
the post-saturation regime and improves the FEL perfor¬ 
mance [13]. 

In order to accomodate diagnostics a realistic undu¬ 
lator design must include periodic break sections, with 
longer breaks adversely affecting performance. This is 
due essentially to three effects. Firstly, diffraction ef¬ 
fects are critical to the performance of a tapered FEL 
particularly for long, multiple Rayleigh length undula- 
tors. While these effects are mitigated by refractive 
guiding inside the undulator, there is no guiding dur¬ 
ing the break sections and the radiation size increases 
exponentially, reducing the field amplitude, causing par¬ 
ticle detrapping and limiting the extraction efficiency. 
Secondly, a break of length Lf, introduces a phase er¬ 
ror Ad/ ~ LbSjA = 2,nXrr} for a particle with relative 
energy offset rj = 5^ with respect to the resonant par¬ 
ticle. Thus longer break sections increase electron phase 
mixing and reduce the bunching factor. Einally as a prac¬ 
tical consideration, for a given total undulator length, 
longer break sections reduce the length of magnetic ele¬ 
ments limiting the electron deceleration and over-all ex¬ 
traction efficiency. To minimize the break length we su¬ 
perimpose the focusing quadrupole field on the helical 
undulator field, similar to the design successfully tested 
in Ref. [14]. One advantage of distributed quadrupole 
focusing is the possibility to operate at small betatron 
beta function, due to the reduced EODO lattice length 


TABLE I: GENESIS Simulation Parameters 


Parameter Name 

Parameter Value 

Beam Energy 

12.975 GeV 

Peak Gurrent 

4000 A 

Normalized Emittances 

0.3/0.3 p m rad 

Average beta function 

5 m 

RMS Energy Spread 

lO”'* 

Bunch Length 

6 fs 

Seed radiation power 

5-25 MW 

Radiation Wavelength 

1.5 A 

Rayleigh Length 

10 m 

Undulator Period 

2 cm 

Undulator Parameter 

3 

Quadrupole Focusing Strength 

26.4 T/m 

Undulator Section Length 

1 m 

Undulator Break Length 

20 cm 

FEL parameter 

1.66 xl0“^ 

3-D Gain Length 

65 cm 



FIG. 2: SASE power and energy spread upstream of the self- 
seeding chicane. The energy spread at 1 GW of power is 3.1 
MeV. 


Lf. This minimizes the transverse beam envelope oscil¬ 
lation = f3avLf/{Pav ~ which also degrades 

the FEL performance m- In our study the undulator 
magnetic field is tapered continuously and the section 
length is chosen to be 1 m, close to the 3-D gain length 
with 20 cm breaks in between. 

Although this kind of undulator has never been con¬ 
structed in the past, the parameters presented in this 
design are similar to what is currently being considered 
for an LCLS-II-like planar superconducting undulator 
with the addition of built in quadrupole focusing |12j . 
A full engineering and tolerance study of this undulator 
is needed before we can be confident that it is a feasible 
option for future high efficiency X-ray FEL facilities. 
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FIG. 3: Bunching factor (left), taper profiles and FEL radiation power evolution (right) obtained from time independent and 
time dependent optimization. The top plots correspond to an input energy spread cr_E,o=3.1 MeV consistent with the SASE 
result from Fig. 1. The bottom plots are an alternate case with (t_b,o=1.5 MeV. In both cases z=0 is after the self-seeding 
monochromator and the input seed power is 5 MW. 


III. SELF SEEDING STAGE 

We consider generating a monochromatic seed pulse 
using self-seeding design based on a single crystal 
monochromator similar to what is currently installed at 
LCLS [TB]. We assume that the tapered section follows 
a self-seeding chicane which completely eliminates the 
beam microbunching, and is preceded by a SASE sec¬ 
tion as shown in Fig.l. The energy spread at the start 
of the tapered section is then determined by the SASE 
FEL process, spontaneous emission losses in the SASE 
section m and the laser heater induced energy spread 
set to suppress the microbunching instability |TH]. We 
perform a simulation of the SASE section assuming an 
initial RMS energy spread due to the linac of 10“^, and 
the results are shown in Fig. 2. 

For effective seeding the radiation power must exceed 
the electron beam shot noise power [T5] : 

^noise ~ ^TTiC /2 ( 2 ) 

by a wide margin. For our parameters this evaluates 
to 36 kW. We therefore decide to start the seeded section 
with 5 MW of power which requires I GW of SASE power 
incident on the monochromator. This sets the length of 
the SASE section Lsase = 13.4 m and the input en¬ 
ergy spread cts =3.1 MeV. As pointed out in Ref. m, 
in practice there is a trade-off between seed power and 


energy spread at the start of the seeded section. Here 
we study the two cases of 5 MW and 25MW seed power 
leaving a detailed analysis of how this trade-off impacts 
the tapered FEL performance for future work. 


IV. TAPERING OPTIMIZATION 
A. Time Independent 

We first obtain the optimal taper profile, maximizing 
the output power for a fixed 100 m undulator length 
in time independent simulations using the three dimen¬ 
sional FEL particle code GENESIS [3D]. The tapering 
law is written as: 

a^{z) = a^o X {1 - c X {z - zo)‘^) (3) 

where the parameters ZQ,c,d are obtained by mutlidi- 
mensional scans which maximise the output power. The 
quadrupole focusing can also be tapered to further in¬ 
crease the extraction efficiency as shown in Ref. [3| but 
that will not be considered in this study. The optimal ta¬ 
per profile obtained from time independent optimization 
is shown in Fig. 3. The tapering order is approximately 
quadratic, which follows qualitatively from the fact that 
in time independent simulations the bunching factor and 
trapping fraction remain nearly constant in the tapered 
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FIG. 4: Maximum radiation power as a function of taper order 
d and taper amplitude ^ for time independent (top) and time 
dependent (bottom) simulations for the case of 1.5 MeV input 
energy spread. 

section and the dominant radiative process is coherent 
emission. The peak output power is 7.3 TW with an 
extraction efficiency of 14 %. It is important to note 
that there is no sign of the taper power saturating in the 
time independent case, which is not the case when time 
dependent effects are included. 


B. Time Dependent Optimization 

Using the optimal taper starting point obtained from 
time independent simulations, we perform time depen¬ 
dent scans over the taper order d and the taper strength 
Attw/aw = c X (Lu, — zqY- As shown in Fig. 3, the 
values of the taper order and taper strength yielding the 
maximum power in time independent simulations are not 
the optimal choice of parameters once time dependent ef¬ 
fects are included. The variation in peak power is more 
sensitive to variations in the taper profile in the time de¬ 
pendent cases. The optimal taper order is weaker than 
quadratic, and is reduced compared to the time indepen¬ 
dent case. This is due to the FEL increased sensitivity to 
particle detrapping when electron beam shot noise and 
multiple frequency effects are included. 

Since the coherent emission power is proportional to 


the product of the number of trapped particles and the 
change in resonant energy (taper strength), a slower ta¬ 
per preserves the trapping for longer, maximising the 
product and the over-all extraction efficiency. In the large 
energy spread case it is important to note that the time 
dependent optimized taper profile has a slower taper or¬ 
der but a larger over-all deceleration rate. This results in 
a worse particle capture in the early stages of the tapered 
section (z=10-50m) but a reduction in detrapping in the 
remainder of the undulator. This can be understood by 
examining the functional form for the resonant phase: 

sinvI/^(z)=x^M (4) 

where x = (2*me*c^/e)(Au,/2As)(l/-\/2[JJ]) is a con¬ 
stant independent of 2 and E(z) is the electric field am¬ 
plitude. The time dependent optimized taper reduces 
|a(„(z)| in the second half of the undulator z=50-I00m, 
maintaining a larger bucket area in the region where the 
amplitude of the sidebands is more appreciable and the 
system is more sensitive to detrapping. From this is clear 
from this that a fully optimized form of the taper pro¬ 
file should have an improved capture rate in the early 
stages with a profile similar to what one obtains from 
time independent optimization, and a slower decrease in 
the undulator field in the later stages when time depen¬ 
dent effects are more appreciable. This requires a more 
elaborate functional form for aw{z) and will be investi¬ 
gated in future work. 

C. Effect of the energy spread 

The input energy spread is a critical parameter for the 
performance of a tapered X-ray FEL. We study this by 
performing the same time dependent optimization for two 
cases both starting with a 5MW seed: the self-seeded case 
with an energy spread of 3.1 MeV and an alternate case 
with 1.5 MeV. In practice the alternate case could be 
achieved by considering a double-bunch system, where 
two closely spaced bunches are separated before the en¬ 
trance to the undulator, the first bunch is sent through 
the undulator to lase producing the seed radiation and 
is discarded prior to the self-seeding chicane. The trail¬ 
ing bunch is transported outside of the undulator and is 
recombined with the seed pulse downstream of the self- 
seeding chicane. In this scheme the seeded bunch would 
have an RMS energy spread set only by the linac and the 
laser heater, around 1.5 MeV for our beam parameters. 

As is evidenced in Fig. 3. the low energy spread case 
achieves a higher peak power, 4.7 TW compared to 3.7 
TW after the time dependent optimization. The taper 
saturation is also delayed due to a decrease in sideband 
induced particle detrapping. In both cases electron emis¬ 
sion into the lower synchrotron sideband mode causes de¬ 
trapping from the high energy region of the stable phase 
space area. Furthermore, scattering of the electrons from 
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interaction with the sideband frequencies causes diffu¬ 
sion and additional particle loss. In the next section we 
discuss theses effects and show how the time dependent 
optimization reduces them. 


V. SIDEBAND INSTABILITY 

The mechanism of sideband generation and amplifica¬ 
tion in free electron lasers can be summarized as follows 
[21]. Firstly, sidebands are generated due to amplitude 
and phase modulations of the electric field, due to the 
trapped particles undergoing synchrotron oscillations as 
they pass through the undulator. Using Maxwell’s equa¬ 
tions in the 1-D slowly varying envelope approximation 
we can write the evolution of the electric field amplitude 
and phase [Hj: 


results in a reduction in particle loss and a delayed ta¬ 
per saturation, both evidenced in the increased bunching 
factor and output power in Fig. 2. In the simple case 
of constant sideband and carrier amplitude the diffusion 
coefficient caused by sideband excitations is proportional 
to the ratio of the power in the sidebands to the power 
in the FEL signal D oc CPs' /Pg with the coefficient C 
depending on the type of sideband spectrum |23j . As 
is also shown in Fig. 4 this is reduced in the time de¬ 
pendent optimized case. The peak power improves by 
1 TW between the time dependent optimized and un¬ 
optimized cases, an overall improvement of 2 %. Despite 
the dedicated time dependent optimization we do not re¬ 
cover the single bucket extraction efficiency unlike results 
previously reported in Ref. pS] . 

A. Effect of the seed power 


/ 

a-s = 7, - 

ZuJsC 

6ks = 


ZujgC a.s 


sin d' 

7 

cosd' 

7 


(5) 

( 6 ) 


where as is the dimensionless vector potential for the 
electric field, uip is the electron beam plasma frequency 
and 4^ is the ponderomotive phase. It is clear from these 
that as the electrons oscillate in the longitudinal phase 
space ('I', 7 ) the gain and the phase shift of the radia¬ 
tion field will be different at different locations in the 
undulator and, due to shot noise in the electron beam, 
at different locations along the bunch. This results in 
a temporal modulation of the radiation field giving rise 
to sidebands displaced from the central wavelength by a 
quantity proportional to the synchrotron period: 


Xs' ~ Ac 


1± 


= A,, 


I± 


I -t- al 


1 / 2 ' 


(7) 


where Lgy is the synchrotron period. Once the side¬ 
bands are generated, the electron oscillations are driven 
by a multiple frequency ponderomotive potential, there¬ 
fore the equations of motion and Maxwell’s equations for 
the electric field, must be modified accordingly. An anal¬ 
ysis of the simplest two frequency model shows that the 
coupled beam-radiation system is unstable and that the 
sideband amplitude will grow from noise for any realis¬ 
tic electron distribution [5T] [55]. When the strength of 
the sidebands exceeds a critical level, electron motion be¬ 
comes chaotic leading to severe particle detrapping and a 
loss of amplification of the FEL signal (53]. Thus, as has 
been discussed by previous authors, suppressing the side¬ 
band instability is the key issue for tapered FEL designs 
[I] , particularly those which are multiple synchrotron pe¬ 
riods in length [53j . 

As is shown in Fig. 4 the time dependent optimized 
taper profile reduces sideband amplitude growth. This 


Increasing the ratio of input seed power to equivalent 
shot noise power in the electron beam minimizes the im¬ 
pact of sideband growth as shown in Fig. 5. We have 
analyzed the impact of this effect numerically by per¬ 
forming the same time dependent optimization described 
above with a 25 MW input seed assuming the same initial 
energy spread of as =1.5 MeV for comparison purposes. 
While time independent tapering optimizations produce 
a peak power of 7.7 TW very similar to the 5MW seed 
case, time depedent optimizations have a much better 
performance, with a final output power of 6.3 TW. This 
is shown in Fig. 6 where we can see the peak radiation 
power still growing after 100 m and the bunching factor 
decaying slowly in the tapered section of the undulator. 
The mitigation of sideband induced detrapping in this 
case is also evidenced by the fractional sideband power 
which remains below 10 % throughout the undultaor (see 
Fig. 4). 


B. Sideband suppression 

In order to reach the single bucket extraction efficiency 
we identify three separate solutions for further suppress¬ 
ing sideband growth which are currently being investigat¬ 
ing. Firstly, as was pointed out originally in Ref. |5T] and 
demonstrated numerically by [53] , increasing the electron 
beam energy spread in the last region of the undulator 
where the sideband amplitude is larger can also reduce 
the sideband growth. The energy spread can be intro¬ 
duced by means of a magnetic delay fine and the interac¬ 
tion of the beam with an external laser inside a short, few 
gain length undulator around the location of exponential 
saturation following the self seeding monochromator. De¬ 
pending on the flexibility of the undulator design one may 
also obtain this additional energy spread by detuning a 
number of undulators around the exponential saturation 
location and allowing the beam to radiate spontaneously 
outside the FEL gain bandwidth HZ]. 
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FIG. 5: (left) Spectrum at half the undulator length and fractional sideband power (right) for the time independent optimized 
case and fully time dependent optimized for 5 MW and 25 MW input seed power and 1.5 MeV energy spread. Sideband power 
grows faster in the time independent optimized case, leading to particle detrapping and early saturation of the tapered FEL. 
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FIG. 6: Time dependent optimized power and taper profile 
for a 25MW seed and 1.5 MeV input energy spread. The 
bunching factor decays slowly due to reduced sensitivity to 
sideband induced detrapping and the radiation power does 
not saturate. 


Furthermore, a wavelength filter with corresponding 
delay line for the electrons could be placed in the tapered 
section in order to select a narrow bandwidth signal [21] . 
This should be done before the sideband power reaches 
the stochasticity threshold beyond which the trapping ef¬ 
ficiency is seriously degraded. Filtering out the sideband 
frequencies will result in a reduced bucket height and 
some initial particle detrapping. This should be com¬ 
pensated by a reduction in sideband induced detrapping 
between the filter and the end of the undulator, thus 
providing a lower limit to the bandwidth of the filter. 
Lastly, we note that the time-independent efficiency can 
be achieved in a fully time dependent simulation by ar¬ 
tificially removing the shot-noise from the electron beam 
distribution. While this is not entirely achievable in prac¬ 
tice, the shot noise suppression method discussed in Ref. 
m could be applied if the scheme can extended to hard 
X-ray wavelengths. 


VI. CONCLUSION 

In this paper we perform the first comparison of time 
independent and time dependent tapering optimization 


for a high efficiency seeded hard X-ray FEL. The compar¬ 
ison is done for an undulator design optimized to achieve 
TW peak powers in the shortest possible distance: he¬ 
lical, superconducting and with built-in focusing. We 
demonstrate that the taper profile yielding the maximum 
power in time independent optimizations does not corre¬ 
spond to the optimal solution when time dependent ef¬ 
fects are included in the simulation. By performing time 
dependent scans of the taper order and the taper strength 
we show that the maximum output power in time depen¬ 
dent mode is achieved with a lower taper order compared 
to the time independent case. The difference is due to the 
increased sensitivity to particle detrapping in the time 
dependent case, mitigated by a slower taper profile in 
the later stages of the undulator where time dependent 
effects are more important. For an input energy spread 
of 3.1 MeV the final output power increases from 2.7 
TW with the time independent taper profile to 3.7 TW 
with the profile obtained from dedicated time dependent 
scans. 

We have also discussed the importance of the trade-off 
between energy spread and seed power at the entrance 
of the tapered undulator section. We show that using a 
“fresh bunch” with input energy spread of 1.5 MeV de¬ 
termined only by the linac we can decrease particle de¬ 
trapping, maintain a larger bunching factor and improve 
the over-all performance. For the same seed power of 5 
MW the maximum output power is 4.7 TW after the ded¬ 
icated time dependent optimization. In a double-bunch 
system the input seed power can be larger without affect¬ 
ing the input energy spread. We have studied an optimal 
case with a 25 MW seed and 1.5 MeV energy spread and 
found that the output power reaches 6.3 TW at the end 
of the undulator, a 12 % efficiency which approaches the 
time independent result of 7.7 TW. 

We identify the sideband instability as the fundamen¬ 
tal time dependent effect which is not taken into account 
in time independent optimizations and limits the extrac¬ 
tion efficiency by causing particle detrapping. Analyzing 
the fraction of energy in the sidebands in the ge =1.5 
MeV case with a 5 MW seed, we show that the fraction 
of energy deposited in the sidebands is below 10 % for 
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70 m in the time dependent optimized taper profile while 
it exceeds 10 % after 40 m in the the time independent 
case reaching 14% towards the end of the undulator. 

While extending the simulation method of Ref. [3] to 
include time dependent effects significantly improves the 
performance of tapered X-FELs the current procedure 
is both time consuming and simulation intensive. The 
form of the taper profile aw(z) needs a more complicated 
functional dependence to optimize trapping in the early 
stages and reduce sideband growth in the later stages 
where time dependent effects play a more important role. 
With the enhanced understanding gained of the critical 
parameters limiting performance, such as the growth of 
the sideband instability, an improved algorithm can be 


developed which acts to directly suppress these effects. 
Such a scheme will be developed in future work. 
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